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Abstract 



We study a correspondence between f(R) model of gravity and a phe- 

c-| ' nomenological kind of dark energy (DE), which is known as QCD ghost 

Qh< dark energy. Since this kind of dark energy is not stable in the context of 

Einsteinian theory of gravity and Brans-Dicke model of gravity, we consider 

two kinds of correspondence between modified gravity and DE. By studding 

the dynamical evolution of model and finding relevant quantities such as, 

equation of state parameter, deceleration parameter, dimensionless density 

5a ' parameter, we show that the model can describe the present Universe and 

also the EoS parameter can cross the phantom divide line without needs to 

any kinetic energy with negative sign. Furthermore, by obtaining the adia- 

batic squared sound speed of the model for different cases of interaction, we 

show that this model is stable. 
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1 Introductions 

It was indicated that our Universe is in a positive accelerating expansion 
phase [1, 2, 3], and the standard model of gravity couldn't explain this 
phenomena. This is a shortcoming for Einsteinian theory of gravity and 
several attempts have been accomplished to justify it. 

One way for explaining this positive accelerating expansion is modified 
gravity. In fact some people have believed this shortcoming of Einsteinian 
theory of gravity is coming from the geometrical part of Hilbert-Einstein 
action and try to modify it by replacing a function of Ricci scalar, f(R), 
instead of R in the Hilbert-Einstein action and so-called f{R) model of 
gravity 4 The various aspect of f(R) models of gravity is investigated in 
[4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20]. 

On the other hand some people looking for the reasons of this accel- 
erating expansion in the matter part of Hilbert-Einstein action. Based on 
this idea, an ambiguous component of matter so-called dark energy (DE) 
is introduced for describing the positive accelerating expansion of the Uni- 
verse. The simplest model of DE is the cosmological constant which is a 
key ingredient in the ACDM model. The ACDM model, which is consistent 
with nearly all observational data, but it faces with the fine tuning problem 
[21]. Also some other DE models have been introducedand are studied in 
[22, 23, 24, 25]. 

Recently, a new kind of dark energy model, so called QCD ghost dark 
energy (QCDGDE) model has been proposed [26]. This model, although 
is un-physical in the usual Minkoski space-tim, express important physi- 
cal effects in dynamical space-time or space-time with non-trivial topology. 
Actually, introducers of this dark energy model, claim the vacuum energy 
of Veneziano ghost flied in effective quantum field theory at low energy, is 
related to density of dark energy which is used for explaining the positive ac- 
celerating expansion of the Universe [27]. They have shown that in a curved 
space time, the ghost field gives rise to a vacuum energy density HAq CD of 
the right magnitude ~ (10~ 3 ev) 4 , where H is the Hubble parameter and A 3 
is QCD mass scale [28, 29, 30, 31, 32]. 

There are two groups which study QCD ghost dark energy model. The 
first group which have introduced it have believed that in the ghost model 
of dark energy, one needs not to introduce any new degree of freedom or 
modify gravity and it is totally embedded in standard model of gravity. 



4 There are another modified gravity formalism which are completely different with 
f(R) model of gravity. 



The second group have accepted the relation of density of QCDGDE as a 
phenomenological model of dark energy and studied the dynamical evolution 
of it in Einstein and Brans-Dicke model of gravity [33, 34, 35, 36, 37, 38, 39, 
40, 41]. The studies of the second group show that the dynamical evolution 
of this model is instable. This means that the QCD ghost dark energy 
model has to be study further. One can study the proposal of the first 
group and consider the gravitational effective field theory at low energy for 
some extra interaction or for some another cutoff for the size of the Universe 
(manifold). And the another possibility for studding this issue is based on 
the second group work with some extra degree of freedom or modify gravity. 
According to second group's view of point, we study QCDGDE model as a 
phenomenological model of dark energy in f(R) model of gravity. 

As was mentioned earlier, modified gravity and DE model of gravity, are 
used for explaining the positive accelerating phase of the Universe and it is 
reasonable there are some correspondence between them. Based on this idea 
we separate some extra terms of modified gravity and define correspondence 
relation between it and a new kind of DE. Moreover, since there are some 
shortcoming in these models when anybody consider one of them alone, 
we studied a combination of modified gravity and dark energy model here. 
Therefore we assume the matter component of the action is consist of cold 
dark matter (CDM) 5 and dark energy. In fact, we combine the extra terms of 
modified gravity with dark energy part of matter component and equivalent 
it with QCD ghost dark energy. Finally, we numerically and analytically 
compute some quantities such as DE density parameter, squared adiabatic 
sound speed, equation of state and deceleration parameter of the model. 

This paper is organized as follows. In Sec. 2, we rewriting field equa- 
tions on the f{R) model of gravity and obtain the equation of motions and 
conservation relation for density energy. In Sec. 3, we make an equivalence 
relation between f(R) model and GDE and obtain the relevant quantities of 
model such as, density parameter, squared adiabatic speed, EOS parameter 
and deceleration parameters. At last we summarize our work and give some 
discussion in Sec. 4. 

2 Description and general properties of the model 

The action of f(R) gravity with general matter is given by 
S = / \J — g d x 



—Z—+L m (ip,g lil/ ) 
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In fact cold dark matter and ordinary matter are denoted by index "m" for simplicity. 



Where f(R) is an arbitrary function of Ricci scalar, R, L m = L^tp^g^) 
is the matter Lagrangian, ip is the matter field, g^ v is the metric of space- 
time and g is the determinant of metric. Here we have assumed 8ttG = 1. 
Variation of (1) with respect of g^ v gives 

iW' - \f9v> + (9^0 - V M V„) /' = T™, (2) 

where prime represents the derivative with respect to the curvature scalar 
R, □ is the covariant d'Alembert operator (□ = V a V a ) and T™ is the 
stress-energy tensor of matter which is defined by 

rpm _ ~2 S{\/—gL m ) 

V" ^ Sgl u, ■ W 

The stress-energy tensor is covariantly conserved, this means that 



V 



R^f - 7;f9^ + {g^n - v„v„) /' 



viz = 0. (4) 



Now we consider a homogeneous and spatially-flat space-time with Friedmann- 
Lemaitre-Robertson- Walker (F1RW) line element 

ds 2 = dt 2 -a 2 {t)(dx 2 + dy 2 + dz 2 ), (5) 

where a(t) is the scale factor. We will assume the Universe is filled out with 
perfect fluids. Then the stress-energy tensor may be given by 

T^ u = (p t + pt) u^u v - g^pt, (6) 

where pt and pt are the energy density and pressure of the fluids and u^ = 
(1,0,0,0) is its normalized four-velocity in co-moving coordinates in which 
u^u^ = 1. So by substituting (5) in to the right hand part of (4) and making 
use of (6), one can arrive at 

Pt + 3H(l+u t )pt = 0, (7) 

where we have used pt = u>tPt- To study the dynamics of DE model in the 
f(R) gravity, we consider a flat Universe with only two energy components; 
cold dark matter (CDM) and dark energy (DE), in which baryon is include 
in the CDM part 6 . This means 

Pt = Pm + PDE, Pt=Pm+ PDE, (8) 



In fact in this work we study a model which consist of f(R) model of gravity and dark 
energy simultaneously. 



so, one can decompose (7) as 



Pm + %Hp m — Q, 

p D E + 3H(l+CJ DE )pDE = -Q- 



(9) 
(10) 



where Q is the direct interaction between two different components of mat- 
ter. One can show that the action (1) is equivalent with [42] 



S 



f(4>)R + $(0) + 2L r 



-g d A x, 



(11) 



where f{<fi) and <&((f>) are the proper functions of a scalar field (f>. So using 
f(R) = f{4>)R + ®((f>), one can rewrite (2) as 



G 



(IV 



L flVl 



where G^, is the Einsteinian tensor and 



T, 



/ IV 



f(4>) 
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and 



rp<p 



-g^H(t>) + (V„V„ - 5^n)/(0) 



It is well known that V^G^ = 0, them from (12) we have V^T) 



//</ 



Since, VT, 



flV 



0, so we can obtain 
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rpm _j_ rp(p 
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The tt component of (15) and Eqs. (9) and (10) gives 

-Q,, 



Pm + %Hp m = Q 

p DE + 3H(p 4 ,+p DE ) = - 
p<p + 3H(p (t> +p (f) ) = 0, 



where 



P<t> = P<p 



f 



3Hf 



[Pm +, 



(12) 

(13) 

(14) 
= 0. 

(15) 

(16) 
(17) 
(18) 

(19) 



As already mentioned in the Introduction, we suppose p^ has the same as 
dark energy role in evolution of the Universe. So we can combine it with 
the dark energy part which is coming from the matter part (the right hand 



side) of Einstein equation, (2). This means that we can rewrite (16 ), (17) 
and (18) as follows 

p m + 3H Pm = Q (20) 

PA + 3H(p A +p A ) = -Q,, (21) 

where 

PA = P<p + Pde, PA=P<t,+ Pde, (22) 

The tt component of the gravitational equations (12), for the metric (5), can 
be simplified to 

"^iffi 1 "^' 1 ' <23) 



where 

Using (23) we arrive at 
H = - 



h(^-3Hf(<p). (24) 



1 



Pt + P<j> + PDE + . 



3 Dynamics of ghost dark energy 



(25) 



In this work we want to consider a correspondence between f{R) model 
of gravity and ghost dark energy. In the ghost model of dark energy, the 
energy density of dark energy is given by pd = aH where a is a constant with 
dimension [energy] 3 roughly of order of Aq CD , where Aqcd ~ 100 MeV is 
QCD mass scale [26]. 7 For this correspondences we have two choices as 
follows 

• We define a correspondence relation as 

Pd = P<t> + Pde, Pd = Pde + P<t>- (26) 

Here p^ and p^ are density energy and pressure of ghost dark energy 
respectively. In this case the relation of conservation for matter and 
dark energy is the same as (20) and (21). 



7 Note that the relevant physical quantities, ujd, q, £l e d, cf, are independent of a. 



• Another possibility for a correspondence relation is 

Pd = P<j> + Pde, Pd = Pde + P<t>- (27) 

In this case the relation of conservation for matter and dark energy is 
as follows 

p m + 3Hp m = Q, (28) 

/ 

Pd + 3H(p d + p d ) = - [p m + p d ] - Q, (29) 

in this case we have two different interactions. Q interaction is a direct 
interaction between CDM and DE and the other term is coming from 
the interaction between matter and geometry. 

3.1 Dynamics of ghost dark energy for pd = Pa 

Some observational data such as, observational of the galaxy cluster Abell 
A586, supports the interaction between DE and CDM [43]. Therefore in this 
section we study the direct interaction between DE and CDM and study the 
dynamical evolution of the model. So in (20) and (21), we assume Q / 0. 
Then by pd = oJdPd we have 

p m + 3Hp m = Q (30) 

p d + 3{l+u d )Hp d = -Q. (31) 

One can rewrite (23), (30) and (31), in terms of dimensionless quantities as 

Q em + Q e d = 1, (32) 

"era H fj - ' 'era T (3 + p)tl\l ern = - — ^ , (3oJ 

2/7 

n ed + lf n ed + (3 + p + 3u D )Hn ed = -^U (34) 

where Q em = p m /3H 2 f (</>), O ec / = pd/3H 2 f((p) are the dimensionless energy 
density of CDM and DE respectively and / = fifH. Using (32) and (33) we 
obtain 

-n ed + — (i - n ed ) + (3 + p)H(i - n em ) = -^- (35) 

and substituting (34) into (35) we have 

2H 

— + (3 + f3)H + 3Hoj d tted = 0. (36) 

ri 



As well, based on the linear relation between p d and H, we have 

Hfl e df = constant, (37) 



Expressing (35) in terms of efolding-number x = In a, and making use (37) 
we have 

" 2 ~ U " ! Q -3+(/?-3)(l-fU. (38) 



-n 



ed~ 



V ed 3fH^ 

Here, the equation of state parameter of the model versus £l ed is 



Ud 



lr3 + 2fiQ- 

"3 [ 2-Sl ed 



0, 



(39) 



where Qq = Q/(3fQ, ed H' i ) and deceleration parameter q = —1 — H/H 2 is 

_ p + 1 - (n Q + 2)n ed 



(2 - n ed ) 

Also we can obtain the squared adiabatic sound speed of our model as 



(40) 



dp d 
dp d 



1 



(3-/3)-(3 + n Q -[3)n ed \ _^_ 

(3 + p) - (3 + n q )n '■"' r 



qj^^ed 



l dn 



ed 



u d . (41) 



For getting better insight we consider the interacting GDE for three different 
forms of Q below. 

3.1.1 Q = 3b 2 Hp d 

In this case Oq = 3b 2 and Eq. (38) reduces to 



n 



ed: ' 



n 



ed 



where 



K 



((3-3-3b 2 )[tt ed + K] 
3-/3 



(42) 



(/? - 3 - 3b 2 ) ' 
and by integrating of it we find 

21n(ft ed ) - (2 + K) Indeed + K\) = (3 - 0)x + C u 

where C\ is the constant of integration and it is given by 

d = 21n(n edo ) -{2 + K) ln(|O edo + K\) 



(43) 



and £l e d is the effective fraction of dark energy at the present time. 
Also Eqs. (39), (40) and (41) reduce to 



Wd 



1 6b 2 - (3 

2 — Q e d 3(2 — £l ed ) 



(3 + 1 (36 2 + 2)Q 



ad 



(2-n fl 



(2-n d 



(44) 
(45) 



n 



ed 



(3 - p) - (3 + 36 2 - /3)O ed 
(3 + /?) - (3 + 36 2 )ft ed j (2 - n ed ) 



A, 



a 



nd) 



-, (46) 



where A\ = (3 + 66 2 — /3)/3. For getting better insight, we compute Ud, 
q and c 2 for /3 = —0.2, b = 0.1 and $7 e ^ = 0.8 and obtain w^ = —1.15, 
q = —0.99 and c 2 = 0.06. It is seen that for this especial choice, the 
equation of state cross the phantom divide line (ujd = — 1) and the adiabatic 
squared sound speed is positive. This means that this model can describe 
the positive accelerating expansion of the Universe and also it is stable. We 
solve numerically (43) and plot it in Fig. 1. Fig. la is the effective fraction of 
dark energy, Q, e d, versus efolding number, x = ln(a), and shows that at the 
early time Vt e d = and the late time is saturated to 1. Also at the present 
time VL e d = 0.8. Fig. lb is the adiabatic squared sound speed versus £l e d- It 
is obviously seen that at the present time namely for VL e d = 0.8, c 2 > 0, and 
then for this kind of interaction the model is stable. 



3.1.2 Q = 3b 2 Hp m 

In this case Q,q = 36 2 (1 — £l e d)/£l e d and equation (38) reduce to 



9. 
where 

By solving (47) we have 



/ ^ ^erf 



ad~ 



Q 



ad 



e(l-fied-l), 



(47) 



£ = 3-/3-36 2 . 



«L 



c 2 e 



& 



|1 — ^ed\ 

where c 2 is the constant of integration and it is given by 



(48) 



C'2 



n 2 , 

edp 



|1 — &ed \ 

9 












' 


- b=0.00 
- - b=0.1 

- b=0.25 




/ " 

/ '■' 

/ /p 
/ " 


- 






1 // 

1 // 

// 
r r t /? t 
// i 

jr i 


" 




- 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 

fi„i 



(a) 



(b) 



Figure 1: (a): This subfigure shows ft ed versus efolding number, x = ln(a). 
(b): This subfigure shows c 2 versus Q ed . We have taken j3 = —0.2, p^ = p\ 
and Q = 3b 2 Hp d . 

where £l e d is the value of effective density parameter of dark energy at the 
present time. And Eqs. (39) and (40) become 



Ud 



e 



26 2 



3(2 - n ed ) vi ed {2 - n 



I'd) 



(3 + 1 - (2 - 3b 2 )n 



(:'(/ 



(2-n«0 

and the adiabatic square sound speed Eq. (41), is 



ai-n ed )[cn 2 ed + Ab 2 (i-n 



ed) 



[(3 + p - 3b 2 ) - (3 - 3b 2 )n ed ] (2 - n ed ) 



(tn ed + 2b 2 



(49) 
(50) 

1 



3^ed(2 — tied, 



-(51) 



We solved numerically Eq. (48) and plot it in Fig. 2. Fig. 2a indicates the 
effective fraction of dark energy, £l ed , versus efolding number, x = ln(a) 
for Q = 3b 2 Hp m . This subfigure shows that the behavior of effective di- 
mensionless parameter of dark energy is the same as Q ed for Q = 3b 2 Hp d . 
Fig. 2b indicate the adiabatic squared sound speed versus Q ed in the case of 
Q = 3b 2 Hp m . It is obviously seen that at the present time namely (? s > 0, 
and then for this kind of interaction the model is stable too. 
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Figure 2: (a): This subfigure shows ft e( i versus efolding number, x = ln(a). 
(b): This subfigure shows c 2 versus Q ec i- We have taken j3 = —0.2, p^ = p\ 
and Q = 3b 2 Hp m . 



3.1.3 Q = 3b 2 Hp t 

In this case Qq = 3b 2 /Cl e d and equation(37) reduces to 



ed~ 



n 



where 



(:(/ 



K 



3 - g - 3b 2 

3-/3 



By solving (50) we obtain 

2 lnfi ed + (K-2) ln(|ft ed - *T|) = (3 - /3)tfx + c, 

where 

c = 21nOo + (K - 2) ln(|fi - K\), 

is the constant of integration. 

l r (3-/?)fi ed + 36 2 1 



3 ^ed(2 — O e d) 



13 + 1 -3b 2 - 2R 



(2 - fid) 
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(52) 



(53) 



(54) 
(55) 
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Figure 3: (a): This subfigure shows ft e( i versus efolding number, x = ln(a). 
(b): This subfigure shows c 2 s versus Q ec i- We have taken j3 = —0.2, p^ = p\ 
and Q = 3b 2 Hp t . 



and 



(3-/3- 3b 2 ) - (3 - /3)ft ed \ (3 - mid ~ 6& 2 (1 - M 



(3 + /3- 3b 2 ) - 30 erf 



■[(3-f3)Q ed + 3b 2 



(2 -a 



(i</j 



3^ed(2 — ^ed) 



(56) 



We solved numerically Eq. (53) and plot it in Fig. 3. Fig. 3a indicates 
the effective fraction of dark energy, £l e d, versus efolding number, x = ln(a) 
for Q = 3b 2 Hpt- This subfigure shows that the behavior of effective dimen- 
sionless parameter of dark energy is the same as Q, e d in the two pervious 
subsection. Fig. 2b indicate the adiabatic squared sound speed versus fl e d 
in the case of Q = 3b 2 H pt- This subfigure shows that at the present time 
namely c 2 > 0, and then for this kind of interaction the model is stable too. 

3.2 Dynamics of Ghost dark energy for pd = p,p + Pde 

In this subsection dark energy is assumed to be without direct interaction 
with the matter part namely Q = . So according to (29) dark energy 
interaction may be with geometry. However the conservation relation of 
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Figure 4: (a): This subfigure shows £l ed versus efolding number, x = m(a) 
(b): This subfigure shows c 2 s versus Q ed . We have taken /3 = —0.2, p d = 
Pde +P<b and Q = 0. 



energy is as (28) and (29), so 

Pm + 3Hp m = 0, 
Pd + 3H (p d + p d ) = j3H [p m + p d ] . 

In this case Eqs. (33) and (34) can be rewrite as 

Iff 

n m + —n m + (3 + p)im m 

ti 

£l ed ~ PHrQ, ed + 3H£l ed + 3u d H£l ed + 2— 9> ed 



H 





(57) 




(58) 


0, 


(59) 


0. 


(60) 



Using Eqs. (32) and (59) one can obtain 



2H 

-n ed + —(i - n ed ) + (3 + p)H(\ - n ed ) = o. 

By combining Eq. 60) with Eq. (61), we have 
2H 



H 



+ 3uj d n ed H + (3 + P)H - (1 + r)pHVL ed = 0. 



(61) 



(62) 



Hence one can achieve in an expression for Q, ed in terms of efolding-number 
x = In a, as 



-n:,(?-J^) + [ i-9..,. 



J e</ 



n 



ed 



3-/3 



0. 



(63) 
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By solving (63) we have 



n 2 



I. o i = ce ^' (64) 

|1 — "ed| 

where £ = 3 — (3 and c is the constant of integration. Eventually, in this 
case, the expression of the EoS, w^, of DE versus £l e( i is specified as follows 



Ud 



and deceleration parameter is as 



1 
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.2 — »'ed 


as 




1 + /3 - 2ft ed 



2-0 



(■</ 



(65) 



(66) 



At last, we can find the squared adiabatic sound speed of the model as 
follows 



1 



n 



ed. 



(3 - /3)(i - n ed ) 
(3 + p) - 3n ed j (2 - n ed ) 



/3(l + r) 



3-/3 

(2-n, 



t.(/j 



(67) 



The numerical results of Eq.(64) and Eq. (67) is plotted in Fig. 4. In this 
case we don't have any direct interaction between CDM and DE. This fig- 
ure shows that the effective dimensionless energy density of dark energy is 
started from at the early time and increased to asymptotic value, 1 at late 
time. Fig. lb indicates the adiabatic squared sound speed versus Q e d and 
shows that at the present time, Q e d = 0.8, c 2 s is positive. So as was expected, 
the model is stable in this case too. 

4 Conclusion 

We investigated two kinds of equivalence between f(R) model of gravity 
and an phenomenological kind of dark energy so- called ghost dark energy 
whose energy density is proportional to Hubble parameter. We studied the 
interacting and non interacting case of model in a flat FLRW background. 
For getting better results, we consider three kinds of interaction between 
matter and dark energy. The obtained results for equation of state and 
deceleration parameter, show that the model can describe the accelerating 
expansion phase of the Universe and also uod can cross the line ojd = — 1 from 
quintessence to phantom model for all cases of model with/without inter- 
action. Also by obtaining the numerical results for effective dimensionless 
energy density of dark energy, £l e d, we find that Q e d is started from at the 
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early time and saturated to its asymptotic value, 1, at late time. We further 
studied the dynamical evolution of the model by considering the adiabatic 
squared sound speed. This quantity is obtained for all cases of interaction 
and without interaction. Our investigation show that the adiabatic squared 
sound speed can be positive and then the model can be stable by a suitable 
choice of parameters. 

At last, we conclude that the dynamical behavior of the model for two 
different correspondence relations, also for all interactions/non interaction 
cases have the same behavior. This model is stable and can describe the 
present Universe. 
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